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We have proposed a new stochastic interpretation of the sudiffusion de- 
scribed by the Sharma-Mittal entropy formalism which generates a nonlin- 
ear subdiffusion equation with natural order derivatives. We have shown 
that the solution to the diffusion equation generated by Gauss entropy 
(which is the particular case of Sharma-Mittal entropy) is the same as the 
solution of the Fokker-Planck (FP) equation generated by the Langevin 
generalised equation where the 'long memory effect' is taken into account. 
The external noise which pertubates the subdiffusion coefficient (occuring 
in the solution of FP equation) according to the formula D a —> D a /u where 
u is a random variable described by the Gamma distribution, provides us 
with solutions of equations obtained from Sharma-Mittal entropy. We have 
also shown that the parameters q and r occuring in Sharma-Mittal entropy 
are controlled by the parameters a and (u), respectively. 

PACS numbers: 05.07.-a, 05.40.-a, 05.60.-k 

1. Introduction 

Subdiffusion is a process where the random walk of a particle is strongly 
hindered by the complex structure of the system. Over the last decades the 
anomalous diffusion process has been observed in many physical systems (see 
[H El [3] and the references cited therein) . One of the most used anomalous 
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diffusion models is the Continuous Time Random Walk model [2] which 
provides a linear anomalous diffusion equation with fractional derivatives. 
Within this model subdiffusion occurs when the random walker waits an 
anomalously long time to take its next step (the mean waiting time is infinity 
and the length of jumps has finite moments). 

The Continuous Time Random Walk formalism provides the relation 

((x-x ) 2 ) = 2D a t a , (1) 

where ((x — xo) 2 ) denotes the mean-square displacement of a random walker, 
xo is its initial location, D a is the anomalous diffusion coefficient measured 
in the units m 2 /s a and a is the anomalous diffusion parameter, < a < 1 
for subdiffusion, a > 1 for super diffusion, and a = 1 for normal diffusion. 
The relation ([T]) has often been used as a definition of subdiffusion. How- 
ever, normal diffusion is described by the stochastic models in which the 
process is assumed to be a Markovian one. As is shown in [1], there is a 
non-Mar kovian process where the relation ([1]) is fulfilled for a = 1. This 
process combines subdiffusion and superdiffusion features, which causes the 
effects of both processes to neutralize each other mutually and provides a 
relation typical of normal diffusion. This important result shows that we 
cannot find a stochastic interpretation on the basis of the relation ([I]). Fre- 
quently, the opinion is expressed that the definition of anomalous diffusion 
should be based upon the stochastic interpretation of this process. 

Let us note that in the system where the relation (pQ) is valid, other 
functions f(t) occur ensuring the relation f(t) ~ t a and which have the 
macroscopic interpretation and are experimentally measured, as, for exam- 
ple the time evolution of the near-membrane layer thickness [5], the time 
evolution of the reaction front in the subdiffusive system with chemical re- 
actions [6] or the functions characterizing subdiffusive impedance [7]. Let us 
also note that there are models which do not have the stochastic interpreta- 
tion (or such an interpretation has not been found yet) but in these models 
we also can find the important characteristics of the system, which satisfy 
the relation f(t) ~ t a . The examples are the models based on nonextensive 
entropy formalism. We do not exclude these models but we are searching for 
their stochastic interpretation. The aim of our paper is to find the stochas- 
tic interpretation of anomalous diffusion model based on the nonextesive 
Sharma-Mittal entropy formalism. 

The description of diffusion can be given by extensive or nonextensive 
entropy (see for example [8j [9] and references cited therein) . As is shown 
in [8] , in the case of the description of diffusion as a / 1 one can only use 
nonextensive entropies. Within the framework of these models, the differen- 
tial nonlinear diffusion equations with derivatives of the natural order were 
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obtained [HI O [TUl 021 HU Q3]. The problem is that such models do not 
have fully satisfactory stochastic interpretation. Until now, attempts have 
mainly been made using the modified Langevin equation for the descrip- 
tion of a diffusion process which is simultaneously described by a nonlinear 
equations [H [21 [151 EH E] ■ The result is that a random force occuring in 
the Langevin equation depends on the solution to the nonlinear equation. 
This situation can be interpreted as an existence of feedback between the 
system and random force which — at least in our opinion — is not easy to 
interpret. We will show that the assumption of such feedback is not neces- 
sary in order to obtain the stochastic interpretation of the diffusion model 
based on nonextensive entropies. 

In our paper we will use a two-parameter nonextensive Sharma-Mittal 
entropy for the description of the subdiffusion process. This entropy has 
a general character; many other kinds of entropies (amongst others Tsallis 
entropy and Gauss one) can be treated as particular cases of Sharma-Mittal 
entropy [18] . First of all we will show that the anomalous diffusion equation 
generated by nonextensive Gauss entropy can be obtained from the Fokker- 
Planck equation obtained from the Langevin equation describing the 'long 
memory' process in which a particle movement is generated by the internal 
noise. On this point we will use the model presented in |19|, 120] . In the next 
step we will assume that the external noise can disturb the subdiffusion 
parameter occuring in the solution of Fokker-Planck equation. We will 
show that the distribution of this parameter connected with the Gamma 
distribution changes the equation and its solution in such a way that we can 
use nonextensive Sharma-Mittal entropy for the description of this process. 
To pursue this goal we will generalise the ideas presented in [21] in the case 
of Sharma-Mittal entropy and the case of subdiffusion. In the next part of 
our work we will only consider a subdiffusion process, i.e. we will assume 
that < a < 1. 



The most general form of entropy we are interested in is the two-parameters 
Sharma-Mittal entropy defined as [8] 



where P is a probability density function of finding a random walker at point 
x at time t. From Sharma-Mittal entropy one can obtain other entropies, 
e.g. Tsallis entropy for q = r, Gauss entropy for r — > 1~ and Boltzmann- 
Gibss-Shannon for r — > 1~ and q — > 1~. For two statistically independent 
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systems A and B entropy satisfies the following equation 

S S m(A + B) = Ssm(A) + Ssm(B) + (1 - q)SsM(A)S S M(B) . (3) 

Thus, for q ^ 1 one deals with nonextensive entropy. 

Entropy ([2]) provides the following diffusion equation El [22] 



dP S M(x,t) 

at 



Q 



SM 



MAr-i) d 2 P r SM (x,t) 
dx 2 



(4) 



where Qsm is the fluctuation strength. The solution to Eq. (JH) for the 
initial condition 

P(x,0) = 5(x-x ) , (5) 

reads [H [22] 



P S M(x,t;x ) = D S Ai(t) 
where {z} + = max{z,0}, 

D SM {t) = 



1 



CsM(t) 



(r - l)(x - x y 



2r(l + g)Q5Af^r,gk r | 2 t 



1 

1 + 9 



C SM {t) = 2{z r D SM {t)f , 

r(r/(r-l)) r>1 

1 , 

1/3 < r < 1 



r-1 r((3r-l)/(2(r-l))) ' 

7T , r 

— r((l+r)/2(l-r)) 
1-r r(l/(l-r)) ' 
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q — r 
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The function 
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provides the relation 

(x - x y 



, r = 1 . 

2 1 

3r - 1 CW(t) ' 



^(aj - x ) 2 ) = j (x- xq) 2 Psm(x, t; x )dx 
Comparing ([1]) with ([6|) we get 
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Q5M = 4rK r>2/a _ M W) ■ (13) 

Thus, the fundamental solution ([6]) fulfils the relation ([1]) only if its form is 
as follows 



P S M(x,t;x ) 



_ (r - l)(x - x f 
2D a (3r - l)t a 



i 

r-1 



^2D a {Zr - l)t a \z r 

(14) 

Let us note that the parameter q is controlled only by the subdiffusion pa- 
rameter a (see Eq. (|12p ). whereas Psm is controlled by three parameters 
a, D a and z r . Let us also note that the parameter r is not related to the 
subdiffusion parameters a and D a . This situation is different from the Con- 
tinous Time Random Walk formalism where the subdiffusion parameters a 
and D a fully determine the process. Let us note, for subdiffusion q > 1. In 
the following considerations we will assume that 1/3 < r < 1 which causes 
functions (|14jl to have an infinite support. 

For r — > 1~ the entropy takes the form of Gauss entropy 

1 _ P («-i) / ^InPdx 
= ^ > (15) 

and Eq. (|14|) is transformed into the following function 

1 (^-^o) 2 
P G (a:,t;xo)= ; =e 4 ^' a , (16) 

which satisfied the following equation 

dP G (x,t) x ^Pg{x^) f 

—df- = aDat dx* ■ (17) 

3. Stochastic approach 

Subdiffusion is a process with 'long memory' which can be described by 
the generalized Langevin equation [191 EQ] 

M dv^ + M j* ^ _ t , W) = R{t) ^ (lg) 

where M is a mass of a diffusing particle, v — its velocity, 7 is the friction 
coefficient, R — the random force with the following correlation function 

(R(t)R(t')) = F t- a . (19) 
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It was assumed that diffusion is caused by the internal noise, thus the 
fluctuation-dissipation theorem can be applied in the form [191 120] 



k,B is the Boltzmann's constant. As is shown in [19] . the equations (flap . (|T9j) 
and (|20p provide the Fokker-Planck equation for a system without external 
forces 

^p(M) _ &p{x,t) 

~^f- ~ ~M Km ~ m] ^^~ ' (21) 

where the functions K and -ff are defined by their Laplace transform L [/(£)] = 
Kp) = I O °f(t)exp(-pt)dt 

K(p) = -— r^r , H{p) = = \, , . (22) 

P + 7(p) P 2 + P7(p) 

Using the formula L\t v \ = T(l + u)/p 1+u , within the limit of small p, 
which corresponds to the limit of long time (in practice the condition t S> 
[MUbT / {FqT{1 — a)T{l + a))} l ^ l ~°^ is enough), after calculations we obtain 

dP(x,t) {k B Tf a ^P(x,t) 

8t F T{l-a)T(a) 8x 2 ' 1 ' 



Equation (|23p has the form of Eq. ()17p . Comparing these equations we find 
the relation between correlation function coefficient Fq and subdiffusion 
coefficient 



D a T(l - a)r(l + a) 



(24) 



Thus, the Langevin equation (fTHj) provides the subdiffusion equation the 
same as the one obtained from nonextensive Gauss entropy formalism if 
the random force fulfils relation ()19p . where Fq is given by (|24p . Then, 
the solution of Eq. (|23p with the initial condition ([5]) is identical with the 
function (fl~6|) . 

Now, we will show that the external noise can disturb the system de- 
scribed by Eq. (|23|) in such a way that this system will be described by the 
solution of the equation obtained form Sharma-Mittal entropy ()14p . To find 
the relation between (fl~4"|) and the solution of Eq. (|23|) we use the integral 
formula 

/ C a-1 e" K dC = P" a r(a) , (25) 
Jo 
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a > 0. Combining Eqs. (fTi|) and ([25]) . after calculations we get 

1 (a-z n ) 2 

= e 4* a OW«)du 



20rt Q (D Q /«) 

= -PsmOM;^o) , (26) 



where r] = j3 — h, f3 = . The function in the square brackets is identfied 
as a Gamma distribution function 

f(u;k,e) = -^-e k u k - l e- eu , (27) 

r(fc) 

with 

fe = p - I = -1±J- , e = /3 - - = 3r ~\ . (28) 
H 2 2(1 -r) ' P 2 2(1 -r) V ; 

Equation (|26p can be rewritten in the following form 

i (a;— xq)^ 

/(«; /3 - 1/2, /3 - 3/2) e ~s^) rf n = P SM (x, t; s ) • 

2\/iTt a {D (X /u) 

(29) 

Equation (|29|) can be interpreted as follows. The subdiffusion coefficient 
is subjected to the external noise which changes the effective subdiffusion 
coefficient from D a to D a /u, where u is the random variable of the Gamma 
probability distribution (|27|) with the parameters given by (|28|). Since it is 
assumed that 1/3 < r < 1, we have j3 > 3/2, k > 1 and O > and we 
also note that f(u, k, 0) — > when -u — )• 0. Thus, large perturbations of the 
subdiffusion coefficient D a have a small probability of occuring. 
The mean value and variation of (1271) reads 



k 1 + r 1 , . < 9 \ A; 1 + r 

<"> = e " ■ << A '"> 2 > = W = (5=1? ■ ( 3 ») 



These parameters fulfil the relation k = + 1, so the Gamma distribution 
is here controlled effectively by one of these parameters fc or 6. It is easy 
to see that 

k = j^-r , © = ^—r , (31) 
(u) - 1 (u) - 1 

thus, the Gamma distribution is controlled here by its mean value alone. 
From Eqs. (J2ZJ) and ([28]) we find 

r = , \ x 1 . 32 
3 (u) - 1 K ' 
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Therefore, parameter r, unsteady up until now, is controlled by the mean 
value of Gamma distibution describing the external noise. Then, Eq. (|14p 
can be rewritten utilizing the subdiffusion parameters a, D a and the pa- 
rameters describing the distribution of external noise, e.g. in the form 



PsM(x,t;x ) 



1 r(fe + l/2) 



^2D a OTTt a T(k) 



(x - x ) 2 
2QD n t a 



-| fc+l/2 



(33) 



where and k are controlled by the mean value (it) 



4. Final remarks 

In our paper we proposed a new stochatic interpretation of the sud- 
iffusion described by the Sharma-Mittal entropy formalism. The other 
entropies (e.g. Tsallis and Gauss entropies) frequently used in modelling 
anomalous diffusion can be treated as particular cases of the Sharma-Mittal 
one. Sharma-Mittal entropy generates a nonlinear diffusion equation with 
the natural order derivatives. We have shown that the solution to the dif- 
fusion equation generated by Gauss entropy is the same as the solution of 
the Langevin generalised equation where the 'long memory effect' is taken 
into account (Eqs. (|18p and (|19p ). The external noise which pertubates the 
subdiffusion coefficient according to the formula D a — > D a /u where u is 
a random variable described by the Gamma distribution, provides us with 
solutions of equations obtained from Sharma-Mittal entropy (f33|) . We have 
also shown that the parameters q and r occuring in Sharma-Mittal entropy 
are controlled by the parameters a and (it), respectively. 

Our consideration has concerned the case of 1/3 < r < 1. Lets us 
note that this case does not contain Tsallis entropy, since Tsallis entropy 
describes subdiffusion for q = r > 1. 
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